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SUMMARY 


around spac^ ^ a P °" erful P um P circulating fluid 

motion of the journal general * „ J ^ , the bearin 8- A small whirling 
channel. The hypothesis that the f luid°f low^i 638 Variation Progressing around the 
the pumped f luid velocity ifgreatfr the L \ , Whlrl Whenever the mean of 
variation wave is discussed and mmn a • -u pe ^ lpheral speed of the thickness- 
bearing instability. « is shomTEat f T Slmple ^P^-^ions of journal- 

hypothesis predicts instability'^onset Correct ly C f or 'unload bearl " 6S the 

overpredicts the onset speed as SrLS™£LSer 0aded h^ 1 " 58 ^ 


INTRODUCTION 


phenomen„n f in^hydrodynamicall^lubricated ^ ,"1“ *h. oil whip 

in 1925 by Newkirk alHarfor (ref If^o ^ 0 ^ f 8S ‘, U ” 3S flrst 
a simple explanation of ST! lighhy leaded W^I TT* eXperlmentS a " d « ave 

= -2 ° f 

When compare^with^conventiona^dynamic stabilitfa^^ysi^otr 0111116 ^ *“ 19 ®°* ' 

c^r f or=e^il b -^ 

wSrSSSi If*™ 3) 

bearing y ls dlscussed ’ first for an unloaded bearing, and then for a loaded 


SOMMERFELD'S ANALYSIS 


The idealized case treated by Sommerfeld (ref 3j in 1904 i 0 

fluid lubricant with uniform viscosity u The fluid fi - 1 d f lth 3n lncoin Pressible 

dimensional; i.e axial fi™, t acosit y T he fluid flow is taken to be two- 

, i.e., axial flow is assumed to be negligible. In addition it is assumed 
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Figure 1* - Equilibrium configuration of journal rotating at rate ft under 
load W. Journal center J displaced distance e with respect 
to bearing center B. 

that cavitation does not occur. The width of the bearing normal to the plane of 
the figure is b. 

Sommer f eld’s principal result is that in the equilibrium position under a 
vertical load W the journal is eccentrically displaced by a distance e in the 
horizontal direction. Because of the eccentricity e the film thickness h varies 
with position 0 around the annular space. For h G << R the approximate relation is 

h = h + e cos 0 = h (1 + a cos 0) (1) 

o o 

where a = e/h Q is the eccentricity ratio. Under the assumptions of laminar viscous 
flow with no pressure variation across the thickness of the film, the only possible 
flow patterns that satisfy the requirements of fluid mechanics in a uniform channel 
of thickness h are combinations of the two basic patterns shown in Fig. 2. Here x 
is distance along the channel and y is distance across the channel with 0 < y < h. 
The fluid velocity (in the x-direction at position y) is denoted by u. The volume 

flow rate across a section of the channel of width b normal to the plane of the 

figure is denoted by Q and 9p/9x is the pressure gradient along the channel. For 
the linear profile at the left of Fig. 2 the pattern depends on the parameter U 
which is the velocity of the upper channel wall (the velocity of the lower channel 
wall is taken to be zero). For the parabolic profile at the right of Fig. 2 the 
pattern depends on the parameter A which is a velocity whose magnitude is four 
times the peak velocity in the profile or six times the average velocity. While the 
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Figure 2. - Component laminar flows In narrow channel with uniform pressure 
across channel. Either linear or parabolic velocity profiles. 


relations in Fig. 2 are strictly correct only for steady state flows with constant 
values of the parameters h, U, and A, the Reynold's theory of lubrication extends 
them to apply to slow variations, both in time t and space x, of these parameters. 

In application to the bearing of Fig. 1, the component flows of Fig. 2 are 
superposed, with x = R0 and U = Rfl, with h given by (1), and the undetermined 
parameter A to be fixed by the requirements of continuity [3Q/36 = 0] and uniqueness of 
pressure [ p ( 0 ) - p(0 + 2 tt)]. The value of A so determined is 


A^ = Rft 


6(1 - a 2 ) 

2 + a 2 1 + 


_L 

a cos 0 


and the total volume flow rate is 


. Rnbh, (3) 

2 + a 

The pressure gradients of the component flows in Fig. 2 are superposed and integrated 
to obtain the pressure distribution p(0) acting on the journal. The resultant of the 
pressures acting on width b of the journal is a force, acting vertically upwards 
through the journal center J, of magnitude 


W = 12TT l J t>R' 


(2 + a 2 )(l - a 2 )* 
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The viscous shearing stresses acting on the boundary of the journal also produce a 
resultant force on the journal but its magnitude is smaller than the pressure 
resultant (4) by a factor of order h 0 /R, and thus may be neglected. The force (4) 
must then be equal and opposite to the applied load W for the journal to be in 
equilibrium. 

The remarkable property of the Sommerfeld bearing model is that the equilibrium 
displacement is at right angles to the applied load. This characteristic implies 
that an unloaded bearing is unstable with respect to slow forward whirling motions 
of the journal. To see this, imagine that an external agent moves the center of 
the journal J of Fig. 1 in a circular path of radius e in the counterclockwise sense 
about the bearing center B. When the journal passes through the position shown in 
Fig. 1 if the motion is slow enough the lubricant flow and pressure distribution will 
be very nearly the same as for the equilibrium configuration shown there, which means 
that the fluid will be exerting a force very nearly equal and opposite to W on the 
journal. This force, in the same direction as the velocity of the journal center J 
in its circular path does positive work on the whirling motion. The Sommerfeld 
bearing thus promotes whirling instability for very slow forward whirling rates. To 
consider more rapid whirling rates it is necessary to extend the Sommerfeld analysis 
to include motion of the journal center J. 


WHIRLING STABILITY OF UNLOADED BEARING 


According to (4) the equilibrium position for an unloaded bearing (W = 0) has 
zero eccentricity (a = 0, e = 0). In this configuration the parameter A^ of 
Eq. (2) vanishes and the lubricant flow pattern is simply the linear profile shown 
on the left of Fig. 2. The stability of this configuration is discussed, first by 
describing Newkirk and Taylor's explanation of half-speed whirl (ref. 1) and the 



Figure 3. - Traveling wave of lubricant film thickness variation circling 
channel when center of rotating journal J whirls at frequency £0 around 
circle of radius e. 
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k^at^cs'ar^di'splajel H 

s rrrj d T“£r sel 

t css advance around ^“^'S.'tL. 


h(0,t) 


h + e cos (0 - cot) 


(5) 


Note that the dependence of h on space 0 and time t is ch a c n f o 

^££,*5 r ular chaaael “ dth - — .1 : STffS 

^irstsi ssrss^. 

(5). The continuity requirement applied to a differential arc of length Rd0 is 


J2.+ b 

R30 + b 9t 


0 


(6) 


s %t^t:%rd lov retalns the iiaear -° fiie ° f 


£ 9h 

2 90 


+ £-° 


(7) 


“esfLYunf™ <5> 18 Substlt “^ ia O) th. result Is that continuity cannot be 


“-f 


( 8 ) 


Newkirlc and Taylor took this to provide analytical verification of the half- 
frequency whirl phenomenon which they observed in a vertical shaft running in a 
bearing with plentiful oil supply. It also appeared to explain the "Oil Resonance" 
peak m response when the rotation rate was twice the natural frequency of the s^tem 

wtere the whlrSne ?r«r 8 1688 s “ lstactor y ln explaining the oil-whip phenonenon 
where the whirling frequency remains at the natural frequency as the rotation rate 

fi is increased, although it was noted that the onset of whipping always “cureTat 
speeds equal to or greater than twice the natural frequency 

th’iSSS^ 

at the bearing to u = R[] at the iournal This fi C 7 varies linearly from u - 0 

f1 „ ® n journal. This flow can be decomposed into a mean 

enco ^ l0C f y P large 1 vo r t ic i ty^ ^ if ifass^ef thaf ^lemlTn till Ifavaf ablfto 
"7 Si invof if gT^esVv afa'tL^^ 

Jf U thf™ C f ref - 2) postulates that energy will b l punpedTw t^ whirl 

around thoV^^ '££ TroZZtlL 

Neutral ZZTZyZZs iZZllZZTZZZ Iq^ UTZl ' 
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hypothesis thus "explains" half-frequency whirls and oil whip for a system which has 
an unloaded bearing with plentiful oil supply. If the system provides little 
restraint on the journal its motion will be primarily determined by the fluid film 
forces acting on it. When a low frequency whirl (w ^ ft/2) is accidently started 
energy will be pumped into the whirl, accelerating the whirl until w gets 
sufficiently close to ft/2 that the energy pumped into the whirl just balances the 
system energy losses in a steady "half— frequency" whirl. If the system provides 
considerable restraint on the journal and only permits appreciable whirling motion 
at a natural frequency the fluid film in the bearing will remove energy from 
accidental whirls at a natural frequency w whenever ft < 2(0. However, when ft > 2<o 
the fluid film will pump energy into any such whirl. Whether or not oil whip 
occurs depends on whether the energy supplied by the fluid film is sufficient to 
overcome the system losses. 

Both of the previous explanations are incomplete in the sense that they do not 
make use of all the physical requirements involved. Both arguments use the flow 
pattern of the underlying steady centered flow and both use the kinematics of a 
small whirling perturbation. The Newkirk and Taylor argument also makes explicit 
use of the continuity requirement for the perturbation but neither explanation 
explicitly invokes any consideration of the pressures developed in the oil film. 

Historically, the first complete dynamic analysis of the whirling stability of 
an unloaded bearing was given by Robertson in 1933 (ref. 4). The development which 
follows is essentially just a linearized version of Robertson's analysis. We 
consider the whirl defined by the eccentricity e (e « h D ) and the frequency u> to be 
a small perturbation on the underlying centered rotation in which the flow 
distribution is simply 


u = Rft y/h , 0 < y < h Q 


(9) 


When the journal is whirling the flow pattern will be a superposition of the two 
profiles of Fig. 2. To fit the conditions of Fig. 3 we take 

u = Rft y/h + A (y/h - y 2 /h 2 ) (1°) 

where h(0,t) is given by (5) and A is to be determined from the continuity require- 
ment (6) and the requirement of single-valued pressure [p(6) = p(0 + 2ir) ] . Using 
a linear perturbation analysis we neglect terms of order (e/h D ) in comparison 
with terms of order £/h 0 and find 

A = -6R ($ - u>) cos(0 - wt) OD 

h Z 


The corresponding total volume flow rate for the lubricant film is 


Q = Rftbh/2 - £R(| - <d)b cos(0 - 0)t) 

= Rftbh / 2 + £ Rbw cos (0 - wt) 
o 


( 12 ) 


The pressure distribution is obtained by integrating the pressure gradient given in 
Fig. 2. To first order in £/h 0 


279 



(13) 


p(0,t) - p o = 12yR 2 ^3 (| - o))sin(0 - wt ) 

h 

o 

journal°displacement 'e'and'hasbthe ^ " 


F = 127Tpeb —r- (~ — Oj) 

h J 1 
o 


( 14 ) 


velocltyby^of “the journal" *1 * “'J h3S th * Sama as instantaneous 

o„ the (j?e d : u r P ^ e ao^ t r^ a ^f z* fluid fiin d ° -* 


FVj = 12Tryb e 2 w(§ - w ) 

h 1 

o 


(15) 


™ S pressure 1S ln h the ?^L t «l“ P rt tUde V* the CO,,,ponent fl °» “1th parabolic profile, 

noJZTiL^t th ~ 

“ hirlS are hlscouraged. The whirling ^re^nc/of ' ‘ 


previously J* c0 ” pa “ d “ Uh the two simplified arguments considered 

previously. The assumption in the Newkirk and Taylor argument that the velocitv 
profile remains linear is equivalent to assuming that the parameter A in Fig 2 
vanishes which according to ( 11 ) implies that (a must equal SI/ 2. Furthermore 
the vanishing of A implies an absence of pressure gradient and consequently an 

fluid C m re , S ^ tln8 f ° r ^ e 30 that at the Particular whirl frequency oo = Si/ 2 the 
iuid film neither retards or advances the whirl. This is the neutral stability 
condition. The heuristic hypothesis (based simply on the flow pattern of the 

^r 8 Cen ^ red /° nfiguration and the kinematics of the whirl perturbation) 
that the mean flow drives the whirl whenever the mean fluid velocity is greater 

than the phase velocity of the whirl is essentially a qualitative statement 
equivalent to the quantitative statement represented by eqn. (15) it happens 

° r a ;\ “"loaded bearing the frequency of neutral Zhl?l is gLn co^Sy by 
the heuristic hypothesis. For loaded bearings it is difficult to see how the 
Newkirk and Taylor argument can be extended to predict any other frequency of 
neutral whirl than oj - ft/2. The heuristic hypothesis is however easily extended. 
It no longer predicts the exact frequency of neutral whirl but it provides useful 
approximations for moderate loads. 


WHIRLING STABILITY OF LOADED BEARINGS 


, . , . th J s section a linear perturbation analysis is made to determine the 
. i n g stability of the equilibrium configuration of Fig. 1. The heuristic 
ypothesis is then applied to predict an approximate value of the whirling 
requency for neutral stability. We consider a small whirling motion of the 

^i: 1 w?:r tered i abOUt 1 the equmbrium position of Fig. I where the journal center 
1 h an 2 ular velocity u m a counter-clockwise sense about a circle of radius 
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e centered on the equilibrium position which lies a distance e to the right of the 
bearing center. In the equilibrium configuration the film thickness h, the parameter 
A]_, the volume flow rate and the load are given by eqns. (1), (2), (3), and 

(4) respectively. When the whirling motion is added the film thickness becomes 

h(0,t) = h + e cos 0 + e cos (0 - cot) (16) 

o 

With this value of h the velocity profile is taken to have the form of eqn. (10) with 
the parameter A to be determined anew from the continuity requirements (6) and the 
requirements of single-valued pressure [ p ( 0 ) = p(0 + 2tt)]. Using a linear 
perturbation analysis with respect to e/h Q (but including terms of all order with 
respect to a = e/h Q ) we find 


A - 


£ 6R 

h 1 + a cos 0 
o 


[fi f (0,t) + a) g (0, t) ] 


(17) 


where 


f (6,t) 


2 

1 - a cos (6 - 0)t ) 

0 . 2 1 + a cos 0 

2 + a 


6a 


(2 + a 2 ) 2 


COS GO 1 


(18) 


g(6,t) 


cos(0 - oot) 


3a 


2 + a 


cos o:t 


Note that (17) reduces to (2) when £ = 0 and to (11) when a = e/h Q = 0. The 
corresponding total volume flow rate is 


Q = 


Q 1 


cRb 


60 — — coswt - wg(9,t) 

(2 + a ) 


(19) 


which reduces to (3) when e = 0 and to (12) when a = 0. The corresponding pressure 
gradient from Fig. 2 when integrated gives the pressure distribution in the 
lubricant film. The resultant force acting on the journal due to these pressures 
has horizontal and vertical components given by 


H - 


12TrpbR 

h 3 (l-a 2 ) 3/2 

o 


1 - a 

2 + a" 


0 - co 


£ sin uit 


( 20 ) 


V = W H 12'n'ybR . 

1 h 3 (l-a 2 ) 322 

o 


1 - a 2 n 0 2 a 2 0 + 

0 - 0 ) + 3a — 


2 + a 


(2+a 2 )cj 1 

(2 + a 2 ) 2 J 


£ COS COt 


which reduce to (4) and (14) respectively when £ = 0 and when a = 0. The forces (20) 
may be decomposed into three forces: the steady-state load W]_ = W, a force 

proportional to [ (1 - a 2 )S7/(2 + a 2 ) - co] which whirls in phase with the velocity of 
the journal center, and an oscillating vertical force whose magnitude is small when 
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a is small . 
cycle; i.e., 


The total work done on the journal by these forces 
the energy per cycle AE imparted to the whirl is 


during one whirling 


AE 


12iT 2 hbR 3 £ 2 

(2+a 2 )(l-a 2 ) 3/2 


4 - 2 a 2 + a 4 

2 + a 2 


ft 



( 21 ) 


The lubricant film forces are once more destabilizing for 
for rapid whirls. The neutral stability whirl frequency 


slow whirls and 


stabilizing 


0) 


4 - 2a 


8 + 2a 



2 2 

3 + (1-aV 0 

2 ? “ 

9 - (1-aV 


( 22 ) 


varies from u - fl/ 2 at a = e/h c = 0 for an unloaded bearing to oo - fl/3 when the load 

W^accordW Jo f22W ^ eCCentr ^ c f t y ratio a approaches unity. The variation of 
according to (22) is represented by the curve labeled A in Fig. 4. 



Figure 4 . - 
ecc en tr i c i 
accordi ng 


Frequency 0) of neutral 
ty ratio: A, accordin 
to heuristic proposal 


stability whirl as function of 
g to complete stability analysis 
(25 1 . 


steady 

( 22 ); 


state 

B, 


w " Urn± "f ao V° the heuristic hypothesis that the lubricant flow pumps energy 
nto the whiri whenever the mean flow velocity is greater than the phase velocity Rw 
of the thickness variation wave, we recall that in the equilibrium configuration 
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of Fig. 1 the film thickness h is given by (1) and the volume flow rate is given 
by (3). The mean velocity 


- ^1 = Rft 1 - a 2 

bh 1 + a cos 0 0 , 2 
Z + a 


(23) 


varies with location around the annulus. To obtain a single velocity to compare 
with the phase velocity of the thickness variation wave, we propose to take the 
root-mean-square average of U. The heuristic criterion for neutral stability then 
becomes 


R 03 


= U 

rms 


Evaluating the integral, we obtain 



co 


n V/ 4 
(1-a ) 


2 + a 




(24) 


(25) 


as the neutral stability whirl frequency according to the proposed heuristic 
hypothesis. The variation of this frequency with eccentricity ratio a is 
represented by the curve B in Fig. 4. The heuristic prediction (25) is not 
identical to (22) obtained from the dynamic stability analysis, but it provides a 
useful approximation for values of the steady state eccentricity ratio that are less 
than a = 0.5. 
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